Quantum non-demolition measurement of microwave photons using engineered 

quadratic interactions 



O 

as: 
< 

OS 



Chunqing Deng, 1 J. M. Gambctta, 1 and A. Lupasjcu 1 

1 Institute for Quantum Computing and the Department of Physics and Astronomy, 
University of Waterloo, Waterloo, ON, Canada N2L 5G7 
(Dated: August 20, 2010) 

We present a quantum electrical circuit with Josephson junctions formed of two anharmonic 
oscillators coupled with an interaction 57? 7I where 71 and 72 are position-like coordinates. This 
type of coupling allows the quantum non-demolition measurement of the energy of one oscillator 
by monitoring the frequency of the second oscillator. Despite the fundamental tradeoff between the 
coupling strength g and maximum photon storage capacity of the oscillators, it is possible to achieve 
high fidelity detection of up to 10 photons over time scale of the order of microseconds. We discuss 
the possibility of observing quantum jumps in the number of photons and related applications. 
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The quantum measurement of harmonic oscillators is 
a generic problem in quantum mechanics. One context 
in which this problem was analyzed is the detection of 
gravitational waves [1]. Another domain of application 
is quantum optics, where the measurement of the number 
of photons in single or multiple modes of electromagnetic 
fields is an essential tool for quantum field state tomogra- 
phy [2, 3], linear quantum computation [4], and quantum 
communication [5] . Most measurement schemes are ham- 
pered by either limited efficiency or large backaction. An 
ideal measurement of an oscillator in the basis of Fock 
states can be achieved by employing a quantum non- 
demolition (QND) strategy [6]. In contrast to non-ideal 
measurements, QND schemes offers interesting prospects 
for preparation of Fock states on demand and application 
to detection of weak forces [7]. 

QND detection of Fock states requires an interaction 
between the observed oscillator and the employed detec- 
tor which depends on the position and/or momentum of 
the former in a non-linear way. Finding a suitable non- 
linear interaction is a difficult task. Up to date, this 
has only been achieved in a few experiments. Exam- 
ples include the measurement of an electron in a Pen- 
ning trap [8], the measurement of microwave photons 
in a Fabry-Perot superconducting resonator [9], and the 
measurement of microwave photons in an on-chip super- 
conducting cavity [10]. Furthermore, other QND mea- 
surement schemes based on nonlinear interactions have 
been proposed for nanomcchanical systems [11] and for 
superconducting resonators [12]. 

In this letter we propose a method for the QND mea- 
surement of Fock states of a simple lumped electrical os- 
cillator formed by a direct-current superconducting quan- 
tum interference device (DC-SQUID) and a capacitor. 
The Fock states are photon number states stored in this 
microwave resonator. We show that it is possible to de- 
sign an interaction between the observed resonator and 
a second resonator which is quadratic in position-like 
coordinates. Combined with a proper choice of reso- 



nant frequencies of the two resonators (i.e. sufficiently 
large detuning) this interaction can be used to measure 
the number of photons in the resonator of interest by 
monitoring the frequency of the second resonator. This 
method for QND detection of microwave photons is dif- 
ferent to the experiment in [10], which was based on using 
the dispersive regime of the Jaynes Cummings Hamilto- 
nian [13], and the proposal in [12], which uses Josephson 
junctions to couple different modes of a single resonator. 
Our method has the advantage of very strong achievable 
coupling resulting in a large signal to noise ratio (SNR) 
for detection of up to 10 photons in the resonator. We 
perform a detailed analysis of the optimization of the 
non-linear coupling used for measurement and the trade- 
off between the strength of coupling and anharmonicity 
of each circuit. 

The electrical resonator that we consider has the sim- 
ple structure shown in Fig. la. It can be seen as an 
electrical LC circuit formed of the capacitance C and a 
SQUID acting as an effective inductor. The SQUID is a 
two-terminal superconducting device formed of a super- 
conducting loop interrupted by two Josephson junctions. 
A SQUID can be used as a sensitive magnetometer, based 
on the property that its critical current depends on the 
magnetic flux $ enclosed in its ring with a period equal 
to the flux quantum $0 = h/2e [14]. For a current in 
the SQUID below this critical current, the SQUID be- 
haves as an inductor, with an inductance that depends 
on the magnetic flux. Therefore, the circuit in Fig. la 
can be seen as an LC resonator, with a resonance fre- 
quency dependent on the magnetic flux enclosed in the 
ring. As such, this system was used for the readout of a 
flux qubit [15]. 

At non-integer values of the applied magnetic flux, a 
finite screening current / s (circulating in the ring of the 
SQUID) arises. This is illustrated in Fig. lb for dif- 
ferent values of the applied flux / = $/$o and ratio 
fi = L/Lj of the geometric inductance of the SQUID 
ring, L, to the Josephson inductance of each junction 
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FIG. 1. (a) DC-SQUID based resonator: the cross symbols 
represent Josephson junctions, (b) (Color online) Screening 
current I s versus transport current I t for / = 0, /3 = (con- 
tinuous blue line), / = 0.4, p = (dashed red line), and 
/ = 0.4, p — 0.2 (dotted black line ), see text for details. 

Lj, with Lj = tfo/I c {fo = ^o/2tt)- For a symmetric 
SQUID this screening current is an even function of the 
transport current I t , and in particular it has a strong 
quadratic component. Therefore, different energy eigen- 
states in the SQUID, characterized by different transport 
current quantum fluctuations, will result in different av- 
erage values of the screening current. A second SQUID 
based LC-resonator, acting as a flux detector, can then 
be used to measure the screening current and therefore 
implement a measurement of the energy, or equivalently 
of the photon Fock states, of the first SQUID. This is the 
basis of the principle of detection. In the following we 
present a full quantum mechanical calculation necessary 
to validate and understand the limits of this qualitative 
argument. 

Fig. 2 shows the full circuit details used to model the 
measured (N) and detection (D) resonators. The two 
SQUIDS are symmetric, with junctions critical currents 
I c ,cn intrinsic junction capacitances Cj jCt , inductance L a 
in each arm (not shown), and magnetic flux f a applied 
to each loop using an external source (not shown), with 
a = N,D. Each SQUID is shunted by a capacitor C a . a 
with a much larger capacitance than the junction capaci- 
tors. The two SQUIDS are coupled through their mutual 
inductance M. A bias current source source I ,N{t) is 
used to control the state of the observed resonator. The 
detection resonator is coupled using a coupling capacitor 
C g to a transmission line of characteristic impedance Zq. 
A wave of amplitude Vj and frequency ujd is sent to the 
detector resonator; the reflected wave amplitude, V a , is 
measured and the result is used to infer the state of the 
N resonator. 

We applied the rules of circuit quantization [16] to de- 
rive the quantum model of the complete circuit shown 
in Fig. 2. The two resonators can be described by six 
degrees of freedom: the gauge invariant phases ji tCt over 
the Josephson junctions (i = 1,2 labels the two Joseph- 
son junctions in each SQUID) and the phases "f a ,a (de- 



FIG. 2. Diagram of the measurement scheme (see text for 
details) . The presence of n photons in the N resonator causes 
an n-dependent shift of the frequency of the D resonator, 
which is monitored in a reflection measurement. 

fined as the flux of circuit theory divided by the reduced 
flux quantum ipo = <I>o/27r [16]) over the shunt capaci- 
tors C a ,a- We use [16] to derive the Hamiltonian of the 
complete system: 

U = ^ T a + ^ Uj t a+UL+H ct T,N+H meaSt z)- (1) 
a=N,D ct=N.D 

T a &U&pC Jta ) + \ P lJ{vl2Cj, a ) + 

\v\ a/(<PoCa,a) is the sum of the electric ener- 
gies in the capacitors in resonator a with p s ,a.i 
Pt, a , and p a ,a the momenta canonically conju- 
gate to the variables j StOI = (71, Q — 72,q)/2, 
7*,q = (7i,a+72,a)/2, and7 ajCe respectively; C a ,N = C a .N 
and C a< n = C a ,D + C g . The Josephson energy in each 
SQUID is given by Uj <a = —2tp I CiC , cos j t ,a cos j a>a . 
The total magnetic energy stored in the two cir- 
cuits is U L = J2a,/3 2( ^0 L a^(7s,a + t/o)(7«,/9 + 

with the inductance matrices L = ( L m Ln) ano - 
V = ( L «+^« Lo ~f LaD ). La, a is a small unavoidable 
stray inductance associated with C a , a - The N res- 
onator is controlled through i?ctr,iv = —foIb,N(t)^ a ,N, 
by using h,N(t). Finally, the probing of the D res- 
onator is done by coupling it to a transmission line: 

-ffmeas,L> = —GgPa,D Vi/C a ,D Vo ■ 

We choose C a>a to be significantly larger than Cj >a . 
As a consequence, the deg rees of freedom 7s, a and 7t,a 
have a much smaller effective mass and therefore they 
can be eliminated from the problem using the Born- 
Oppcnhcimer approximation. Assuming that these fast 
degrees of freedom remain in the ground state (which 
is justified in view of the experimental parameters we 
choose), we determine the effective potential energy of 
the slow degrees of freedom by minimizing the potential 
energy with respect to j Sta and jt,a- The potential en- 
ergy terms Uj, a and Ul in Eq. (1) are replaced by U' Ja 
and U' L , respectively, which now only depend on "f a ,a- We 
calculate the effective energy in the first order approxi- 
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mation with respect to j3 a costt/q,, with j3 a = L aa /Lj ta 
and we find: 



Uj,a 



2cp Q COS 7T f a COS J a , a 



J .J,a 



(2) 



and 



— (Lj sin7rfcos7 a ) L(L 7 sin7rfcos7 a ) 



2 

■ L 



(Lj 1 cos7rf sin7 a ) T L'(L J 1 cos7rf sin7 a ). (3) 



Lj = (Yl'J and f = ( / N / ° D ) are the Joseph- 
son inductance matrix and bias flux matrix respectively, 
where we have defined the Joscphson inductance Lj >a = 
fo/Ica- 1 a = ( 7" r> ) i s a column vector contain- 
ing variables J a ,a- We define cos7 a = (msvc) and 

/ sin 7„ jv \ 
Sln 7a - (sin 7? ,J- 

After applying the rotating wave approximation, the 
Hamiltonian of the coupled system has the form 

n' cff /h =[lo d + A D {a)a - l)]a f a + [lj n + A N (tfb - l)]b f b 
+ J{a)b + atf) + 2e D (t)(a i + a) + ga^atfb 

(4) 

where a and b are the annihilation operators for the 
D and N resonators respectively, ££>(t) = eocoswoi, 
with ed the driving amplitude, and we have left out 
the H c tr,N in Eq. (1). The resonance frequencies are 
given by oj a = ^p a /Lj ia C a , a , with p a = 2 cos irf a + 
J2p TJJ sin7r/ Q sin tt fp - cos 2 nf a . The other pa- 
rameters in the model are given by 



_ k^/u D uj N ^(3 D p N cos nf D cos nf N 
J — - (o) 



Z^PnPd 
k^/uj D Lj N sJ /3 D f3 N sin tt/d sin -nf N 



(6) 



4n„ 



' (\ cos nf a + -P a sin 2 irf a ) (7) 

,aPa «J ^ 



with /c = M I V LdLn the inductive coupling constant 



(|fc| < 1), and n D 



t RQCa, a /8iT with Rq = h/e 



the quantum resistance. The maximum number of pho- 
tons, n maXjQ , is the number of photons corresponding to 
a energy equal to approximately the modulation height 
of the potential energy for each SQUID; this is the max- 
imum photon storage capacity for each resonator. Equa- 
tion (6) shows that the optimal value of g is at best 
y/oJD^N / y/nmax,Nn-ma,x,D; multiplied by a factor of the 
order of the unity. With given resonant frequencies (fixed 
values of and uin) 9 can be increased at the expense of 
n max .D, however this leaves the fidelity of measurement 
unchanged (see discussion below). An augmentation of g 
is possible by reducing n max ,iv- 

To implement the measurement of the N resonator, 
the D resonator is continuously monitored through the 



coupled transmission line. We introduce the decay rate of 
the resonator k — +Z )} } where Z is the 

characteristic impedance of the transmission line. This 
decay rate is tunable by merely changing the coupling 
capacitance C g . The system is described by the master 
equation 

R = -i[H, R] + k(clRcJ - \a)aR - \Ra ! a), (8) 

where R is the quantum state of the combined detector 
and system. We assume that that we can neglect the sclf- 
anharmonicity of the D resonator and that \u>jj — lun\ S> 
\J\ so that the linear interaction can be ignored (it only 
shifts the frequencies of the resonators). We then follow 
the procedure outlined in Ref. [17]: the reduced state 
p = Tr[i?] = J2mn Pmn\m){n\ of the N resonator is found 
to be p mn = p mn (0) expH Jo H mn dt'-S Q l T™ n dt'] where 



H r 



H™ n + g(m - n)Re[a n a* 



■pmn I _ 



-d t \a r 



(9) 



with H™ n 



(u>n — An)(iti — n) + An(iti 2 — n 2 ) and 
d n = —i&D — i(A + gn)a n — na n /2. Here A = ujd — &d 
is the detuning between the driving and resonance fre- 
quency. In (9) the first equation represents the Hamil- 
tonian evolution of the N resonator which comprises of 
both bare evolution (first term) and the ac-Stark shift in- 
duced by the detection resonator. The second equation 
represents dephasing: it comprises of a term proportional 
to the information extracted by the detector (see below) 
and a second term which can be both positive and neg- 
ative as it represents information flowing between the 
detector and system that has not been extracted by the 
monitoring. 

In order to quantify the detection efficiency we consider 
the problem of distinguishing a n from a n + 1 photon 
state in the N resonator. In this case the measurement 
rate is 



K\a n +\ 



(10) 



This can be interpreted as the leakage rate of the D 
SQUID resonator multiplied by the distance in phase 
space of the two coherent states corresponding to n and 
n + 1 photons in the N resonator. With a detuning 
A = — g(n + 1), the value of T m is optimized for a given 
driving strength amplitude s~d- The driving strength can 
be raised to the point where the number of photons in 
the D resonator reaches the maximum photon number 
n max ,fl. In this optimal condition, the measurement rate 
is 

F% n+1 = 4nmax - Dg2 [1 - exp(-te/2)] 2 . (11) 

K 

For a homodyne measurement, the signal to noise ra- 
tio for distinguishing n from n + 1 photons is related to 
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the measurement rate through the relation [17] SNR = 
(Jo ™ \A?r m (i)dt) 2 /T m , where T m is the integration time 
and r\ is the detection efficiency for the wave reflected off 
the D resonator, in practice limited by the available am- 
plifier. Using Eq. (11), we find 



SNR = r)An maKi D(gT m )' 



\kT„ 



2(1 



-«T m /2\j2 



(«r m ) s 



(12) 



For a given measurement time T m , the maximum value 
of the SNR, SNR max = 0.32 77n max , D (gT m ) 2 , is attained 
when k = 4.3/T m . T m should be chosen to be of the 
order of the relaxation time T n+ i iU of the N resonator, 
as further increasing will corrupt the signal to noise 
ratio. The optimization of the measurement amounts to 
maximizing <? 2 n maX! 7j, based on Eq. (6). The result of 
the optimization is only weakly dependent on the res- 
onance frequencies, as long as linear coupling is effec- 
tively suppressed (\ujd — <^n\ 3> | J\). High-efficiency mi- 
crowave transistor based amplifiers with i] w 1/20 are 
available over a wide frequency range in the GHz do- 
main, which allows for flexibility in choosing ujjj given 
any a; at. Note: with the advent of quantum limited am- 
plifiers [18], it will be possible to have 77 < 1. To illustrate 
the optimization procedure, we choose W£>/27r = 7 GHz 
u>n/2tt = 4 GHz. Even though Eq. (6) indicates that 
when fo and /jv tend towards 1/2 the value of g can 
reach its upper limit, we choose Jd — 0.45 and /jv = 0.45 
in order to reduce the effects due to the asymmetry of 
the junctions and also to avoid too small values for the 
maximum numbers of photons. We also assume the mag- 
netic coupling constant to be k = 0.7. We then proceed 
to optimizing PdPn /n maXt NPNPD, under the constraint 
2(3 a cos 7r f a -C 1 (we take 2/3 a cosirf a to be less than 
0.1 and fixed values of lojj and u>n). Note: pn and pn 
are implicit functions of /3d, /3/v, Wd, and ujn. In the 
optimization procedure we choose to consider values of 
the ratio Lrj/L^ between 0.5 and 2. The reason for 
this is that very dissimilar values of the loop inductances 
make it difficult to reach large values of the magnetic cou- 
pling constant k with simple circuits. The optimization 
of the SNR over /3d, Pn, an d Ld/Ln, leaves still one 
free scaling parameter, which we take to be Ljn. Wc 
set Lj t N to be 1.1 nH, corresponding to a critical current 
7 c ,jv = 300 11A. 

In Fig. 3a we show the optimal SNR, for 77 = 1/20, 
versus the Lrj/L^ ratio. With a value of T m = lp,s, 
consistent with relaxation times of superconducting res- 
onator with Josephson junctions [19], the SNR is above 
2000, which results in a fidelity for distinguishing n = 
from n = 1 above 99% [20]. The SNR for distinguishing 
?i from n + 1 scales as 1 /n 2 due to the 1/n dependence of 
the relaxation time T Il+ i „. Therefore we estimate that 
the observation of up to 10 photons in the resonator can 
be done with an SNR larger than 20 (state dependent) 
which should result in fidelities larger than 87%. Fig. 3b 




FIG. 3. (Color online) (a) The SNR as a function of L D /L N . 
(b) The maximum number of photons that can be store in 
each resonator n max ,jv and n m ax,D versus Ld/Ln- 



shows that 7i max jv is large: up to 30 photons can be 
held in the N resonator, and 7i maXj £> is large enough for 
the linear approximation to hold. We note that to reach 
the values of the SNR plotted in Fig. 3a, the parame- 
ters g/2ir, Au/2w, Am/2tt, and J/27T are of the order of 
10 MHz. At this value of J, linear coupling is negligible. 

In conclusion, we presented a new method for quantum 
non-demolition detection of photons in a microwave res- 
onator, based on a optimized non-linear interaction. We 
find that large signal to noise ratio can be obtained for 
up to 10 photons, assuming resonator relaxation times of 
the order of microseconds. The quantum non-demolition 
nature of the detection enables the monitoring of quan- 
tum jumps and related applications to detection of weak 
electromagnetic fields [7] . 

This work was supported by an NSERC Discovery 
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